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ABSTRACT 

New aspects of the slow solar wind turbulent heating and acceleration are investigated. A 
physical meaning of the lower boundary of the Alfven wave turbulent spectra in the solar at- 
mosphere and the solar wind is studied and the significance of this natural parameter is demon- 
strated. Via an analytical and quantitative treatment of the problem we show that a truncation 
of the wave spectra from the lower frequency side, which is a consequence of the solar magnetic 
field structure and its cyclic changes, results in a significant reduction of the heat production 
and acceleration rates. An appropriate analysis is presented regarding the link of the considered 
problem with existing observational data and slow solar wind initiation scenarios. 

Subject headings: Physical data and processes: Magnetohydrodynamics (MHD), Physical data and pro- 
cesses: Waves, Physical data and processes: Turbulence, Sun: solar wind, Sun: magnetic topology. Sun: 
atmosphere 



1. Introduction and the physical content 
of the problem 

An understanding of the heating and acceler- 
ation processes of the solar atmosphere expand- 
ing into interplanetary space as a continuous so- 
lar wind and solar eruptions such as coronal mass 
ejections is one of main targets of solar physics. 
It has been realized that the observationally de- 
termined characteristic Reynolds numbers provide 
suitable conditions to develop MHD turbulence in 
the strongly magnetized solar (stellar) wind. The 
interplanetary space represents a 'natural labora- 
tory' able to sustain turbulent flows and allows 
one to test results obtained from basic physical 
groun ds. In the latter context, for instance, re- 
cently IChen et aL I (I2OIOI) studied an anisotropy of 
the solar wind turbulence between the ion and 



dence in favor of the Alfven-cyclotron dissipation 
heating process based on the obs ervations of he- 
lium a nd hydrogen temperatures; ISahraoui et al 



( 2009t ) reported the first direct determination of 
the dissipation range of MHD turbulence in the 
solar wind at the electron scales. These and 
many other similar studies naturally include and 
are complemented b y astrophysical consequences: 
Marino et al. I (I2OO8I) gave an estimation of the rate 



of turbulent energy transfer, which can contribute 
to the in situ heating of the wind; iLuo fc 
( 2010l ) revealed the scale-dependent anisotropy of 
the wave power spectrum. In addition, suprather- 
mal ions and plasma wave spectra upstream of in- 
terplanetary shocks driven by coro nal mass ejec- 
tion e vents have been analyzed by iBamert et al 
(I2OQ8I) . 



electron scales. iKasper et al.l (|2008[ ) provided evi 
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Already in the 70s of the last century the turbu- 
lent cascade of Alfven waves has been suggested 
as an efficient process to transfer energy from 
large-scale waves to small-scale ones, which are 
more easily subject to damping by different trans- 
port processes possibly operating in the plasma 
flow. Fundamental work, particularly for the so- 
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lar context, has been done by Tu an d collaborators 
(|Tu et alJIloillTu fc Marschlll996f ) by developing 
basic aspects of the theory of MHD turbulence 
in the solar wind. More recent works demon- 
strate a consistency of the wave heating process 
with the observed ten iperature and velocity espe- 
cially in coronal holes ( Hu et al]ll999l : Vainio et al 



20031 ) ■ We know that si milar processes can operate 



in stellar atmosphe r es (Narain fc Agarwa l 199 



Ulmschneider et all I2OOII: lElfimov et alj [200 



Therefore, proper modeling of the turbulence is 
of broad astrophysical interest. 

When the modeling concerns the equatorial re- 
gions (even at solar minimum characterized by the 
'double structure' of slow and fast wind) the wave 
heating scenario is rather complex because of sev- 
eral reasons. Firstly, at low latitudes the mag- 
netic field topology drastically differs from that of 
the polar regions and, secondly, there is clear ob- 
servational evidence for multiple sources sustain- 
ing the observed he ating and accelerat ion rates 
in equatorial regions (jKasper et al.ll2Q07l ) through 
the solar cycle. Many analytical and numerical 
models use artificial ad-hoc functions (see , e.g. . 



van der Hoist et all l2007t iKleimann et al.l UQM) 



mimicking the observed latitudinal variation of the 
solar wind structure at solar minimum. The mod- 
eling is even more complex close to solar maximum 
when the solar wind becomes a complicated mix- 
ture of 'channels' of fast (above coronal holes) and 
slow (above coronal streamers) wind regions. 

It is well known that the solar (and plausi- 
bly many stellar) atmosphere(s) is (are) abun- 
dantly populated with Alfven waves (s ee, e.g., 
ICirtain et aDl2007l : IdePontieu et al.ll2007l ). 

A stochastic nature of the wave excitation and 
propagation processes naturally led in the past 
to the notion of MHD turbulence in the con- 
text of the heating and acceleration of the solar 
wind. It is now conventionally accepted that an 
engagement of the clearly observed Alfven wave 
spectrum into the turbulent cascade s represents 



one of the important pro cesses (e.g.. IChen et al 



20101: iMarino et all l2008t lOughton fc Matthaeus 



20051 ) for a transmission of the wave energy from 



large (non-dissipati ve) to small (damping) scales 



(jKasper et al.ll2008l) 



No n-modal cascading was suggested bv lShergelashvili 
( 20061 ) as an alternative process that can be specif- 
ically relevant for the boundaries of coronal holes 



et 



in the regions where the transition from fast to 
slow wind occurs. Actual realizations of such 
processes and their proper testing in comparison 
with observational data is one of the paramount 
challenges of modern solar (and generally stellar) 
physics. 

In order to perform a proper analysis of the 
solar wind formation processes one should model 
the complex and very dynamic structure of the 
solar atmosphere and outflowing wind. There 
are many processes with a wide range of char- 
acteristic spatio-temporal scales which govern the 
complex topological configuration of the magnetic 
field. The aspects playing a dominant role in 
the distribution and properties of the Alfven wave 
sources are the time-varying streamer belt and 
the related heliospheric current sheet. At so- 
lar minimum the streamer-like structures and the 
current sheet are concentrated in the vicinity of 
the ecliptic plane. This configuration reproduces 
the double (slow-fast) structure of the solar wind. 
Namely, at low latitudes the slow wind is gener- 
ated (with characteristic velocity at solar mini- 
mum of about 400 km/s), while almost all of the 
rest of the solar surface in both hemispheres is cov- 
ered by huge 'polar' coronal holes giving origin to 
the fast wind (of about 800 km/s). There is a rel- 
atively sharp boundary between these two zones. 
At this stage there are no spatial constraints on 
the wave sources as the central streamer structures 
are rather wide and maintain a larger range of fre- 
quencies within the streamers. In general, how- 
ever, we suppose that there is still some restric- 
tion on the frequency domain (as we model below) 
within the streamer belt even at the cycle mini- 
mum manifested by the presence of the slow wind 
at low latitudes. However, with the advance of the 
solar cycle this topological picture changes drasti- 
cally. In particular, the latitudinal width of the 
streamer generating zone widens gradually reach- 
ing very high latitudes at solar maximum (up to 
70° and even further). This zone becomes more 
densely populated with the streamer-like struc- 
tures with significantly smaller transverse spatial 
scales compared to those at solar minimum. As 
a result, characteristic wavelengths (frequencies) 
of waves excited in these structures can be bound 
from above (below). Observations clearly evidence 
tfiat t he shape of the solar w ind evolves accord- 



ingly (iRichardson et al.l 120011 ). The velocity pro 
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file becomes 'homogeneous' (characteristic speed 
600 km/s) throughout the streamer zone and the 
transition to the high speed wind (800 km/s) oc- 
curs only within small polar caps to which the po- 
lar coronal holes are confined. 

The second aspect contributing to the shap- 
ing of the Alfven wave spectra are active regions 
hosting other open field structures, which mainly 
can contribute to the high frequency part of the 
spectrum. The resulting magnetic structures as 
sources of the waves are known as active region 
sources ( Liewer et al 1 12004 ). The proper and ana- 
lytically consistent modeling of this factor requires 
a separate effort bringing the issue beyond the 
scope of the current paper. Here we point out 
at least some implications how the contribution 
from the active region sources could be modeled, 
and a rigorous study of the issue will be pub- 
lished elswhere. In the present paper, which is 
based on the simple analytical framework for the 
solar wind wave tu rbulent heating developed by 
Vainio et al.l (|2003l ). we focus on the effect of the 
frequency truncation from below, which must be 
expected to have an effect on the overall process of 
the wave cascading and damping when this lower 
ferequency boundary is located within the inertial 
range. The range of Alfven wave frequencies ob- 
served in the solar atmosphere is determined by 
wave excitation and damping. The wave damp- 
ing has received great attention in the past and is 
well elaborated in the related literature (including 
the citations above), and the analyses carried out 
in those works clearly indicate the significance of 
the high-frequency (small spatial scale) boundary 
position, which is determined by the interconnec- 
tion between the spatial scales of waves and damp- 
ing processes at work. In order to make clear the 
subject of the analysis that we are presenting in 
this paper we should state that, we focus here on 
the low-frequency boundary of the domain based 
on physical grounds, and we treat it as a natural 
parameter playing a significant role for the rela- 
tive efficiency of turbulent processes in different 
regions of the solar atmosphere. The locations of 
the latter depend on the transversal spatial scales 
of the magnetic field structures and related wave 
excitation sources, which are responsible for the 
creation of the wave spectrum. This issue requires 
more attention because its significance usually has 
been diminished in existing models. This is the 



main goal of the proposed analysis, and therefore 
in what follows we give a mathematical formula- 
tion of the above statements and parametrize the 
modeling in terms of the relevant physical quanti- 
ties using a notation that is conventional for mod- 
els of the wave spectral power diffusion due to the 
turbulent cascade in the solar atmosphere. 

2. The model equations 

A quantitative model of the contribution of the 
MHD turbulence into the total heat balance of the 
expanding solar atmosphere has been developed 
and condensed into the balance equation govern- 
ing a steady-state s patial configura tion of the spec- 



tral wave power P (|Tu et al.lll984r ) 



OF 



(1) 



where u and va are the background and the Alfven 
velocity, respectively, and F represents a 'cascad- 
ing' of the spectral power because of the presence 
of the turbulent cascade. For an ex plicit form of 
the ca scading function F see, e.g., IVainio et al. 
(|2003l ). This equation usually is generalized to 



its dynamical counterpart to ob tain the te mporal 
evolution of the spectral power (|Hu et al.lll999i) : 



f)P P fiF 

— + V-[{u + va)P] + -{V-u)^~— (2) 

Following the conventional notation we define the 
pressure produced by waves as: 



Pn 



1 

8^ 



fo 



Pdf. 



(3) 



Here, fo and fn are the 'boundaries' of the consid- 
ered wave spectrum. These two parameters play 
key roles in the wave heating models and they in- 
clude links between the Alfven wave spectrum fre- 
quency domain and the physical conditions in the 
environment where those waves propagate and in- 
teract. In particular, fo represents a minimal fre- 
quency at which the wave energy is injected into 
the cascade and it contains information on the 
characteristic spatial and temporal scales of the 
wave sources, while fn is a maximum frequency 
beyond which the turbulent spectrum is truncated 
because of strong damping and it manifests prop- 
erties and typical scales at which the considered 
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wave dissipation processes are active. We turn to 
this issue again below. 

Taking the integral of equation ([2]) w.r.t. wave 
frequencies within the mentioned frequency range 
fo < f < fa , after straightforward manipulations 
and taking into account that 

VPdf = V / Pdf] +P(/o)V/o-P(/h)V/h 

ifo yfo J 

(4) 

as well as 

r"l?d/ = F(/H)-F(/o) (5) 

Jfo "J 

one arrives at the familiar equation for the wave 
pressure: 

^ + V.[(^^+^/^K,] + ^(V.^/) + ^ = (6) 

with heating contribution Qw from the Alfven 
wave spectrum, which we write in general form 
as: 



47r An 

(7) 

In the standard models of the solar wind wave 
heating main attention has traditionally been 
given to the spectral truncation limit fn, and 
corresponding explicit expressions (derived on the 
ground of the underlying microscopic transport 
processes in a plasma, say cyclotron damping etc.) 
for this frequency have been given in the related 



literature (see, e.g. iHu et aLlflOoi . Even a gra- 



dient of this limit along the radial distance from 
the Sun has been taken into account manifesting 
itself in the erosion (shrinking from the high fre- 
quency side) of the turbulent spectrum. However, 
a contribution from the lower frequency limit /o 
has always been omitted. It was assumed that the 
characteristic frequency, where the actual wave 
spectrum begins, is constant - implying a van- 
ishing gradient of fo in the expression ([7]) ~ and 
at values significantly below the inertial range of 
the spectrum within which the MHD turbulence 
operates effectively and as a result waves mainly 
follow the WKB behavior. This latter assumption 
removes the wave power cascading term i^(/o) = 
and leads, therefore, to a simplified versi on of ([7]) 
as given, for instance, in Hu et al. ( 19991) . In the 



Pifo 



latter work fo is even set to zero with the corre- 
sponding physical reasoning in the fast solar wind. 

We agree that, in principle, the above argumen- 
tation can be valid, but only for the case when 
the lower boundary of the spectrum always re- 
mains outside the inertial range. However, we ar- 
gue here that this assumption cannot be true for 
the slow solar wind regions populated with curved 
and isolated islands of magnetic structures even 
if both terms containing fo in expression ([7]) re- 
main vanishing because of no cascading to the cut- 
off frequency from lower frequencies (leading to 
F{fo) = and P(/o) = 0). The effect of the spec- 
tral truncation still plays a role when fo appears 
within the inertial range, as it systematically en- 
hances the wave power deficit with increasing dis- 
tance from the Sun. Many theoretical and numeri- 
cal models of the solar wind include physical quan- 
tities set artificially as discontinuous functions of 
latitude for the sake of mimicking the real struc- 
ture of the solar outer atmosphere. However, if one 
wants to develop a consistent model for the evolu- 
tion of the solar wind, then gradual changes with 
heliographic latitude in the topology of the mag- 
-^■^(^jcjipl^^aji^jCp -responding scales of the Alfven 
wave sources should be taken into account. There- 
fore, our main intention is to study wave turbulent 
processes in connection with the variable source 
properties and probabilities of their appearance 
during the solar cycle. In this study we use sim- 
ple modeling of the mentioned variability, in order 
to create a framework for further more complex 
modeling of the dynamical processes in the solar 
atmosphere responsible for the formation of the 
solar wind. Understanding both the methodology 
how to do this and possible ways of realizations 
linked with the actual spatial distribution of the 
observable physical quantities represent the goals 
of the present paper. 

3. Solution of the wave transport equation 

3.1. A class of solutions with a low- 
frequency cut-off 

For the subsequent analysis we require a solu- 
tion of the wave transport equation that includes a 
lower frequency boundary of the wav e power spec- 
trum. Like in the reference model (Vainio et al 
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Fig. 1. — Plot of the function g{x) — I{x, 0) — H{1 — x)/x (solid line). The dashed line shows the part of 
the spectrum that has been truncated. 



20031 ) we start from the equation: 



dr dx 



(8) 



where, the function / and the inde pendent vari- 



ables r and x are defined like in Vainio et al 



(120031 ). i.e.: 



2/3 



(9) 



1/2 



ar . 



y3(r') 



ne{r') 



(10) 

where we have chosen the normalisation frequency 
as fn — fo- The above equation ^ describes a 
simple 'dynamics' of the dimensionless wave power 
I in the (x, T)-space and can be solved by the 
method of characteristics, implying the following 
set of equations: 

dl = 0, (11) 

dr = dx/I^^^. (12) 

As in the reference model, this implies that the 
characteristic "time" variable (we refer to this 
variable as time only formally, it actually cor- 
responds to the heliocentric distance r) satisfies 
t' = T and that the "Lagrangian" differential of 
dl — vanishes so that I = Ic = const, as well as 
X + I^/'^T ~ Xc — const, along stream lines in the 
(x, T)-space. 



Therefore, to find the solution one should look 
for it in the form: 



where 



/ = g(.T + /l/V), 



g(x) =/(x,0) 



(13) 



(14) 



determines the shape of the spectrum at the solar 
surface t — 0. For the sake of clarity we invert the 
analysis and search the solution x{I,t) of EqfT2l 
which is: 

(15) 



x{I,t)=F{I)-I^/^t. 



with F(I) being the inverse function of g{x) and 
T < I^^/'^F{I) determining a characteristic dis- 
tance of the wave transport, where the waves reach 
certain small scales a; — >■ and dissipate. In order 
to make our current derivation directly compara- 
ble with the reference model we define the spec- 
trum at the surface as 

gix) = I{x, 0)^H{l-x) x^ ^H{l-x) go{x) 

(16) 

where H denotes t he Heaviside function, q < 5/3 
( Vainio et al ]|2003[ ). the subscript indicates func- 
tions defined as in the reference model, and the 
value X = 1 corresponds to the frequency fo where 
the spectrum is cut abruptly because of the deficit 
of low-frequency waves in the spectra, incorpo- 
rated by the factor H{1 — x). The function g{x), 
which is shown for g = 1 in Figure [TJ corresponds 
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to the following spectrum at the surface: 



fo 



f 



(17) 

The function F{I) can be represented in this 
case in the following form: 



J po(/,ro). 
^ ' F{I) 



1 if < / < 1 
7 if / > 1 



In the limit x <^ I-^^^t one has: 



1 



if < / < 1 



/— if/>i 



(18) 



1 

7? 



if r > 1 
if T < 1 



(23) 
(24) 

(25) 



Using the solution ([T5|) we arrive at a set of alge- 
braic equations governing the relation between /, 
X and r within the entire domain under consider- 
ation: 

r i_ [^(/,t) + /i/2^] =0, T</-i/2 



- [x (/, t) + = 0, T < if/> 



and 

P{f,r) 
if < / < 1 
1 



Po{f,r)=PwKB 

(26) 



(19) 

It is straightforward to show that for x ^ /^/^t 
the approximate solution is I{x, t) w I(x, 0) and 
that we obtain the same WKB behaviour as in 
the reference model. The only difference is that 
the spectrum P(/, r©) is truncated from below at 
the frequency f = fo- 

Now we turn to the other limit x ^ I^^^t. In 
this case one finds: 



2(39-5) 
T 3(3-9) 



if T > 1 

if r < 1 



(20) 



which can be readily rewritten as: 

Po(/,r)(|^)^=AvKB(/,r)(|)"' 



PiLr) 



Poif,r)=PwKBif,r)(^ify 



(21) 

2 

where fg = foTi-^ is a generalized breakpoint fre- 
quency. 

3.2. Comparison v^rith the reference model 

In the previous subsection we have derived a 
class of solutions of the wave transport equation, 
parametrized by the value of the exponent g, tak- 
ing into account a truncation of the spectrum to- 
wards low frequencies. In order to perform a direct 
comparison with the quantitative analysis carried 
out in the reference model, we briefly list the rel- 
evant equations for the solution for g = 1 (see 
Figure [1]) : 

H{l-x) 



(22) 



where fc = /o/t is the breakpoint frequency as 
defined in the reference model. 

There are several conclusions one can draw 
based on the above derivation. In particular, the 
solutions given in the reference model are inde- 
pendent of the frequency normalization constant 
/„, while the new solution explicitly depends on 
it. As the most convenient value we have chosen 
fn = fo- Further, we can see that the truncation 
of the initial spectrum results in a significant mod- 
ification of the spectra at all heliocentric distances 
and, consequently, one can expect a related modi- 
fication of the heating as well as acceleration rates. 
Exactly the quantitative investigation of the latter 

ts-the subject of the following sections, 
if r > 1 

4. v^i^h\e lower boundary of the fre- 
quency domain 

We start our analysis from the concept of struc- 
tured magnetic fields within the streamer zone as 
opposed to the non-structured 'straight' magnetic 
fields in coronal hole regions. It is apparent that 
the characteristic rate of curvature of the mag- 
netic field should depend on the statistical distri- 
bution of the magnetic structures fed by means 
of two main known contributors - shearing due to 
the random velocity fields and random emergence 
of the magnetic flux. The topological structure 
of the magnetically governed zones, of course, also 
depends on the phase of the solar cycle which leads 
to the time variable strength of the magnetic field 
and to a corresponding variation of the internal 
structure of the streamer structures. 



if r > 1 



if T < 1 
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As we have mentioned above, the presence of 
the lower frequency boundary in the spectra con- 
tains information on the scales and distribution 
of wave sources across the solar disk. This in- 
formation can be imported into the model from 
the observations of the distribution of magnetic 
structures apart from the turbulent cascade pro- 
cess itself. In order to calculate properly the lower 
frequency limit we have to establish a relation be- 
tween the location of the magnetic structures de- 
veloped during the solar cycle and the scales of 
wave sources. 

Before we calculate the physical quantities 
determining characteristic temporal and spatial 
scales of wave sources it ha s to be mentioned that 
(I2OO3I) 



we use 



Vainio et al 



as a reference model 
and throughout the following text we adopt eval- 
uations of physical quantities taken from the sub- 
section 2.1 of that work. 

We formulate a simple model of the observed 
variability of several physical quantities assuming 
the corresponding distribution to be related to the 
wave source probability distribution over latitude 
within both streamer zones and coronal holes. For 
the sake of simplicity we use the smoothed Heav- 
iside function: 



Hix) = 



1 



1 



-2kx 



-,x e (-cx),oo) (27) 



where k is the parameter determining the sharp- 
ness of the transition. The dynamical changes 
in the magnetic field topology (as described in 
the introduction) should impose constraints on the 
available scales Sc (measured in Mm) of the wave 
sources, which we approximate in the following 
manner: 



Sc^ 



300 - 200|sin(f^t') 



(28) 



here 1^2 = (77r/18)| sin(7rt'/ll)| is the location of 
the transition from the zone of finite /o to that of 
a vanishing one. At the solar minimum (t' = 0) 
i?2 = represents the fact that there is no trunca- 
tion of the spectrum even if streamers are at the 
equator. We calculate the lower frequency bound- 
ary of the spectrum (using the Alfven speed) as 
/o — VAifQ)/ Sc to obtain four curves of /o vs. lat- 
itude, each corresponding to four different stages 
of the solar cycle: t' = 0, 1.5, 3.5 and 5.5 years. 
We plot the corresponding curves of /o in Figure [21 



Panel A, for the mentioned stages of the solar cy- 
cle. 

Further we write V = u + va sis a, function of 
time: 

V{t',^) = Vi{t') + {800 ~ Vi{t'))H{d ~ ^i), (29) 

where i9i = 7r/12+(57r/18)| sin(7rt7ll)| is the vari- 
able boundary of the velocity transition, t' is time 
measured in years, k — 50, V and Vi are measured 
in km/s and 



Vi{t') = 400 + 200 



sm 



11 / 



(30) 



Corresponding curves are shown in Figure [21 Panel 
B, in the same order as described above. 

In order to demonstrate how the above scenario 
operates in practice and to prove the validity of 
the approach we perform further a quantitative 
analysis based on the proposed updated model of 
the solar wind turbulent heating. 

5. Quantitative analysis and proof of the 
concept's validity 

We aim to perform an appropriate quantita- 
tive study of the problem in order to obtain a 
consistent proof of the proposed concept. As we 
study the possibility of the presence of lower fre- 
quency cut-offs in the Alfven wave spectra above 
streamer-like structures, it is natural to assume 
that the locations of such restricted sources are 
linked with the latitudinal extent of the streamer 
generating zone and its variation with the solar 
cycle. 

In this paper we show how the frequency cut- 
off affects the wave spectrum as well as the heat- 
ing and acc e lerati on rates. Therefore, like in 
Vainio et al.] (|2003[ ). we use V = 400 km s^^ as 
a target velocity. Using this value we derive four 
characteristic values of the cut-off frequencies: 1) 
/o = 10^® Hz. This is the case recovering the one 
in reference model (case of no truncation) as in this 
case the cut-off frequency remains outside the iner- 
tial range along the entire domain of computation 
so that there is no modification of the spectrum. 2) 
the other valuas are calculated and shown in Panel 
A of Fig. [3] as vertical lines at /o = 1.3 • lO^'' Hz. 
(solar minimum, t' — (solid line), /o = 1.8 ■ 10""^ 
Hz {f = 1.5 (dotted line)), /o = 3.0 • lO'^ Hz. 
(<' = 3.5 (dashcd-dotted hue)) and /o = 4.0 • 10"^ 
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20 40 60 80 20 40 60 80 



■& Heliographic latitude (degrees) d Heliographic latitude (degrees) 

Fig. 2. — Panel A: Sketch of the latitudinal distribution of the lower boundary of the frequency spectrum. 
The curves correspond to four different phases of the solar cycle: t' = 0, 1.5, 3.5 and 5.5 (years counted from 
minimum activity) labeled 1 to 4, respectively. Panel B: Sketch of the observed distribution of the total 
speed V vs. latitude. The four curves are shown in the same order as in Panel A. 




f [Hz] Normalized radial distance 



Fig. 3. — The linestyles correspond to four different stages of the solar cycle in the following order: t' = 
(thin solid line), t' = 1.5 (dotted line), t' = 3.5 (dashed-dotted line) and t' = 5.5 (dashed line). Panel A: 
The scaled wave power spectra (fP) vs. frequency at 0.3 AU. The thick solid line is the complete spectrum 
of the reference model. The vertical lines mark the cut-off frequencies /o at four different stages of the solar 
cycle. Panel B: Ratio of the spectral powers, at the corresponding dissipation frequencies, calculated using 
the considered frequency domain truncation and the one taken from the reference model: The difference 
between the reference values (thick horizontal line ) and the updated turbulence spectra (the curves) is a 
result of their low-frequency truncation within the inertial range. 
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Hz. (solar maximum, t' — 5.5 (dashed line))). In 
the same panel we plot the spectra at 0.3 AU cor- 
responding to these five cases. In panel B of Fig.|3] 
we plot the ratio of the spectral powers calculated 
using the considered frequency domain truncation 
and the one taken from the reference model. The 
points of break at which the curves corresponding 
to the updated spectra start to deviate from the 
reference values (the thick horizontal line) mart 
the distances at which r — 1 (/c = /o). These 
curves were obtained using updated spectra acu- 
rately calculated in this pap er and given by ex - 



pression (case q — 1). In Vainio et al. ( 20031 ) 



it has been shown that there is only a small differ- 
ence between convective and diffusive formulations 
of the flux function F. We concentrate on the first 
approach and consequently use: 



F = 27rC 



V B 



(31) 



The physical quantities in the latter expression are 
defined as follows: 



P{f,r) 



-PwKB (/,?■) 
l + [///c('-)]'/' 



if /o < /c 



(32) 



with 



PwKBif.r) = P(/,rQ)— , (33) 

BqVaq 

representing the part of the spectral power out- 
side the inertial range where waves follow predom- 
inantly a WKB behaviour and where 

P{f,r^) = ep^H(^l-[^l^yy (34) 

and fc{r) is the breakpoint frequency at which the 
Kolmogorov type of the spectrum starts to prevail: 

1 ^^u. r C^r').Ar') 



fc(r) 



re 



y2 



ne{r') 



(35) 

Here ep is a dimensionless spectral parameter de- 
termining the intensity of wave excitation. Its 
value in the case of absence of additional wave 
sources is set to 5 • 10^^. In general, a more com- 
plete model of the solar wind must include also 
angular components of velocities, however for the 
sake of simplicity we assume those components to 



be vanishing. C^(r) = aai{r) is a model parame- 
ter, where a is a cascading constant and 



ai 




Tq < r < lOrg 
r > lOr, 



(36) 







with ao another dimensionless parameter deter- 
mining the ratio of outward and inward propa- 
gating wave intensities. Its value when additional 
sources are absent is set to 0.05. 

We have also performed computations for the 
solar wind accelerati on due to the wave pressure 
gradient Vp„. As in Vainio et al. ( 2003() . we use 
explicit the expression: 



VA Ar. 



VAG 



A 



■PwQ 



2VA 



A{r') 
VA{r') 



Q{r')dr', 
(37) 

where A c>c is a cross sectional area of the 
open field structure. The results of calculations 
are shown in Fig. |4] Panel B. Again, the thick solid 
curve shows the radial distribution of the wave 
pressure gradient —{mpne)~^dpwo/dr correspond- 
ing to the case of no cut-off and, thus, the reference 
model. The calculations for —{mpne)~^dpw/dr 
are done for the above-mentioned four epochs of 
the solar cycle shown by the thin curves with the 
same order of the linestyles as in Panel A. It is 
seen that the acceleration rate decreases signif- 
icantly (according to these calculation by 15-25 
percent) within the first ten solar radii from the so- 
lar surface, while it monotonically decreases in ab- 
solute value beyond this distance. To summarise 
this part of our investigation we conclude that the 
truncation of the available wave spectrum from the 
lower frequency side leads to a significant decrease 
in both heating and acceleration rates and, there- 
fore, this proves the main message of this paper 
that the effect of the lower frequency boundary /o 
should not be abandoned in the slow wind models. 

The analysis above lets us conclude that if only 
one kind of wave sources would operate in the 
equatorial zone then the appearance of more and 
more streamer-like open structures with increas- 
ing solar activity - leading to a gradual decrease 
of their transversal spatial scales - would result 
in significant cooling of the wind plasma above 
the distance of 10 solar radii and deceleration of 
the wind stream even within that zone (what is 
actually observed outsi de the sunspot generatin g 
zone at high latitudes ( Richardson et al.l l200ll )). 
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Howe ver, observations evidence (jKasper et al 
20071 ) that there clearly operate two kinds of 
wave sources with relative significance at differ- 
ent stages of the solar cycle. It is reasonable 
to assume that these growing number of small- 
scale magnetic structures in the vicinity of active 
zones increases the probability of reconnections, 
on the one hand, and leads to creation of addi- 
tional small-scale open field structures anchored in 
active zones, on the other hand. It is convention- 
ally presumed that those reconnection processes 
shake these open structures and, thus, represent 
so-called active region sources of additional wave 
power. A significant number of such small struc- 
tures has a reasonable potential to compensate 
the losses of energy arising from the truncation of 
the lower frequency part of the spectrum. In other 
words, even though for low frequency modes the 
equatorial part of the atmosphere is not accessible, 
increased wave power is nevertheless available in 
the high-frequency part of the spectrum above the 
cut-off frequencies /q. One could ask here: what 
should be an actual value of the higher frequency 
wave power to compensate the losses in thermal 
and mechanical energy production rates obtained 
above? As we mentioned above this issue stays 
outside the scope of the current study. It should 
be noted that observ ed wave spectra (e.g. see 



Tu k Marschlll996[ ) contain a combined wave 



power supplied by the sources from entire coronal 
holes (including polar regions) and some possible 
local sources not located on the Sun, but in the 
solar wind. This is why, possibly, the truncation 
of the power spectra as shown in Figure |3] is not 
observationally resolvable so far. However, large 
error bars in the spectra below the frequencies 
10~^ Hz maybe the indication that the significant 
decline in the wave power for low frequencies could 
be observable with higher spectral resolution. 

6. Discussion and conclusions 

As is conventionally accepted, the physical 
grounds underlying the appearance of the fast 
and slow winds are substantially different. The 
fast wind originates in the unipolar magnetic con- 
figurations of coronal holes with open field line 
topology, while most of the indirect measurements 
of the active region evolution during the solar cy- 
cle, like detections of the helium abundance in the 
solar wind, indicate that multiple processes for 



slow wind heating and acceleration should oper- 
ate in the different proportions at different phases 
of the solar cycle (Kas per et a l. 2007). There is 
the process related to the heliographic current 



sheet and the streamer belt ( Einaudi et al. 2000l 



20011 ). which at the eq uatorial region is n iost ac- 



tive at solar minimum (jKasper et al.ll2007[ ). With 
the uptrend of the activity cycle the significance of 
the streamer belt contribution decreases gradually 
reaching some smaller but finite rates at maxi- 
mum. At the same time more and more streamer- 
like structures appear at high latitudes covering 
practically the entire solar disk at maximum ac- 
tivity with relatively small-scale structures. These 
structures lead to the cut-off frequencies of the 
wave spectra generated within those streamers re- 
sulting in reduced acceleration and heating of the 
plasma outflow. This scenario is in good agree- 
ment with the observations at high latitudes. 

A detailed numerical validation of the presented 
paradigm of the wave spectra formation is beyond 
the primary scope of the current paper where we 
focus on the physical grounds of the underlying the 
concept. While more extended numerical studies 
via direct simulations will be published elsewhere, 
here we demonstrate in general terms a correspon- 
dence of the above-stated scenario with the exist- 
ing data and modern understanding of the solar 
wind acceleration scenarios. 

At the solar minimum the sources located in the 
current sheet and streamer belt dominate and the 
dynamics of the process is govern ed by the am- 



bient low frequency Alfven waves (jEinaudi et al 



20011 ) for which the equatorial region is accessible 



to some extent because of the more regular spa- 
tial structure of the magnetic field below the cusp. 
The structure of the magnetic field resembles a 
dipole field at that stage. With the change of the 
magnetic field topology the streamer zone spreads 
to very high latitudes and, at the same time, new 
open field lines appear which are anchored in the 
active regions and are fed from the latter with the 
modified spectrum of Alfven wa.ve dis turbances 



(jLiewer et all 120041 : iKasper et al.l 120071 ) . There 



fore, the excited wave spectrum, following the spa- 
tial scales of the modified sources, becomes trun- 
cated from the low frequency side. Besides, with 
the increasing number of these new active mag- 
netic structures, the rate of confinement of the 
plasma at low heliospheric distances leading to a 
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Fig. 4. — Panel A: Calculated heating rates Qq (reference model, thick solid line), Q (modified heating using 
the four truncated spectra shown in Fig. [3l Panel B: Calculated wave pressure gradients for the reference 
model — (mpTie)"-^ {dpwo/dr) (thick sohd line) and for the truncated spectra —{nipTie)^^ {dp^/dr). The 
linestyles correspond to Panel A. 
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significant modification of the density profile may 
also grow. 

In fact, the scenario developed in this paper is 
an attempt to formulate a mathematical scheme 
for a modeUng of the slow solar wind through- 
out the entire latitudinal domain and throughout 
the solar activity (;yck;. Of course, our considera- 
tions are based on substantial assumptions regard- 
ing the fine structure of temporal and latitudinal 
variation of physical quantities. More rigorous as- 
trophysical modeling is needed to explore further 
specifics of the Alfven wave frequency domain in 
the solar atmosphere that should comprise, e.g., 
the north-south asymmetry of the sunspot distri- 
bution, deviation from normal latitudinal distri- 
butions, or non-equilibrium components of the ob- 
served statistics. The results can have far reaching 
impact on the understanding of the solar and, in 
general, stellar wind origin. This aspect is also of 
interest in the stellar wind context in connection 
with the currently emerging field of stellar wind 
dynamics and interaction with extrasolar planets. 
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of the European Union Research Training Network 
Solairc (MTRN-CT-2006-035484). Financial sup- 
port through the Research Unit 1048 (projects FI 
706/8-1/2) funded by the Deutsche Forschungs- 
gemeinschaft (DFG) is also acknowledged. The 
work has also been financially suppoerted within 
the framework of the FP7-PEOPLE-2010-IRSES 
project no. 269299 - SOLSPANET. The authors 
are thankful to Jens Kleimann, Stefaan Poedts 
and Carla Jacobs for valuable discussions. We are 
grateful to the referee for very constructive com- 
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